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ABSTRACT 

We present optical spectroscopy of MWC 656 and MWC 148, the proposed optical 
counterparts of the 7-ray sources AGL J2241+4454 and HESS J0632+057, respec- 
tively. The main parameters of the Ha emission line (EW, FWHM and centroid 
velocity) in these stars are modulated on the proposed orbital periods of 60.37 and 
321 days, respectively. These modulations are likely produced by the resonant inter- 
action of the Be discs with compact stars in eccentric orbits. We also present ra- 
dial velocity curves of the optical stars folded on the above periods and obtain the 
first orbital elements of the two 7-ray sources thus confirming their binary nature. 
Our orbital solution support eccentricities e ~ 0.4 and 0.83±0.08 for MWC 656 and 
MWC 148, respectively. Further, our orbital elements imply that the X-ray outbursts 
in HESS J0632+057/MWC 148 are delayed -0.3 orbital phases after periastron pas- 
sage, similarly to the case of LS I +61 303. In addition, the optical photometric light 
curve maxima in AGL J2241+4454/MWC 656 occur —0.25 phases passed periastron, 
similar to what is seen in LS I +61 303. We also find that the orbital eccentricity is 
correlated with orbital period for the known 7-ray binaries. This is explained by the 
fact that small stellar separations are required for the efficient triggering of VHE radi- 
ation. Another correlation between the EW of Ha and orbital period is also observed, 
similarly to the case of Be/X-ray binaries. These correlations are useful to provide 
estimates of the key orbital parameters P rb and e from the Ha line in future Be 7-ray 
binary candidates. 

Key words: stars:emission-line, Be - stars: individual: MWC 656 - stars: individual: 
MWC 148 - gamma rays: stars - X-rays: stars. 



1 INTRODUCTION 

Cherenkov telescopes are opening a new era of discoveries 
with the detection of large populations of 7-ray sources rang- 
ing from galactic objects (supernovae remnants, pulsar-wind 
nebulae, compact binaries) to dista nt blazars and starburst 
galaxies (jHinton fc Hofmann 112003 ). In particular, the class 
of 7-ray binaries has attracted great attention in the last few 
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years. They are characterised by the presence of a compact 
star orbiting a late-O/early B companion. Either if particles 
are accelerated in a microquasar jet or in the collision of 
two winds, the high density of UV photons from the mas- 
sive star provides the necessary environment for the produc- 
tion of GeV-TeV emission through inverse Compton scatter- 
ing. There are currently three confirmed 7-ray binaries with 
High Energy (HE, E>100 MeV) and/or Very High Energy 
(VHE, E>100 GeV) emission modulated on the o rbital pe- 
riod: PSR B1259-63, LS 5039 and LS I +61 303 (|Paredes I 
2011). Three more candidates have been recently proposed 
but their binary nature awaits confirmation. These are 1FGL 
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J1018.6-5856 (|Corbet et al. Il201lt >. AGL J2241+4454 and 
HESS J0632+057. The latter two sources are the subject of 
this paper. 

AGL J2241+4454 is a point-like 7-ray source detected 
by th e AGILE satellite above 100 MeV (|Lucarelli et al. I 
2010). The error circle of the satellite contains the Be 
star MWC 656 (= H P 215227 ) which was proposed by 
I Williams et"aH i|201(t ) as the likely optical counterpart. The 
same authors suggest a spectral classification B3 IVne+sh 
and an orbital period of 60.37±0.04 d based on Hipparcos 
light curve data. By fitting the spectral energy distribution 
in the UV and B— band a distance of 2.6 ±1.0 kpc is derived. 
This, combined with its high Galactic latitude b — —12°, im- 
plies a runaway star, another evidence supporting its binary 
nature. 

On the other hand, MWC 148 (=HD 259440) is 
a B0 Vpe star which has been proposed as the op- 
tical counterpart of t he TeV source HESS J0632±057 
(Ah aroman et al.ll2007l). The claim is based on the detec- 
tion of radio JSkilton et"ai1l2 009) and X-ray emission asso - 
ciated to MWC 148 jHinton et al]|2009l ; lFalcone et alfcoid ) 
with similar spectral indexes and variability as observed 
in PSR B1259-63, LS 5039 and LS I +61 303. Recently, 
extended and variable radio emission at AU scales, show- 
ing similar morphology to these three 7-ray binaries, has 
been reported (Moldon et al. 2011a; see also Moldon et al. 
2011b, Ribo et al. 2008 and Dhawan et al. 2006). Lower 
limits to the orbital period of > 54 d (Fal cone et al. 1 120101 ) 
and >200 d (|Casares et al.ll201ll) were proposed based on 
X-ray and optical emission line variability. Further, the ab- 
sence of significant radial vel ocity shifts in the B0 star sup- 
ports a long per i od > 100 d l|Aragona et al.ll2"010h . Finally, 
iBongiorno et al.l (|201 lh report the presence of strong X-ray 
flares in Swift/XRT data, modulated with a period of 321 ± 
5 d which strongly advocates for its binary nature. 

Alternative scenarios to the binary hypothesis invoke 
an isolated star where high energy emission is produced in 
a wind driven shock. In particular, Bp stars are known to 
posses strong surface magnetic fields which efficiently con- 
fine the s tellar wind and may ac celerate particles up to TeV 
energies l|Townsend et alj|2007h . The best way to test the 
binary model and solve the system geometry is through a 
radial velocity study and hence we have embarked in an op- 
tical spectroscopic campaign of MWC 656 and MWC 148 us- 
ing several telescopes. Our analysis demonstrates that these 
stars are indeed the optical companions of High Mass X- 
ray /7-ray Binaries responsible for the HE/VHE emission of 
AGL J2241±4454 and HESS J0632±057. In the following 
sections we present the observations, our best constraints to 
the orbital parameters and discussion of results. 



2 OBSERVATIONS 
2.1 MWC 656 

We observed MWC 656 using the Fibre-fed RObotic Dual- 
beam Optical Spectrograph (FRODOspec) on the robotic 
2.0m Liverpool telescope (LT) at the Observatorio del Roque 
de Los Muchachos between 23 April and 28 July 2011. The 
spectrograph is fed by a fiber bundle array consisting of 
12x12 lenslets of 0.82" each which is reformatted as a slit. 



The spectrograph was operated in high resolution mode, 
providing a dispersion of 0.35 A pix -1 and spectral resolv- 
ing power R ~ 5500 in the blue arm while 0.80 A pix -1 and 
R ~ 5300 in the red arm. The spectral coverage was 3900- 
5215 A and 5900-8000 A respectively. A total of thirty two 
600s spectra were obtained with the blue arm and sixty four 
290s spectra with the red arm. One blue and two red spectra 
were obtained on each observing night. A log of the obser- 
vations is presented in Table [T] The FRODOspec pipeline 
produces fully extracted and wavelength calibrated spectra 
with rms<0.1 A above 4400 A. The analysis presented in this 
paper has been performed with the FRODOspec pipeline 
products. 



2.2 MWC 148 

MWC 148 was observed with the Intermediate Dispersion 
Spectrograph (IDS) attached to the 2.5 m Isaac Newton 
Telescope (INT) at the Observatorio del Roque de Los 
Muchachos on the nights of 20-23 October 2008. A total of 
35 spectra were obtained in the blue spectral range (3900- 
5500 A, unvignetted) using the R900V grating in combina- 
tion with the 235 mm camera and a 1.2" slit to provide a 
spectral resolution of 80 km s _1 (FWHM). Nightly aver- 
ages were produced from the individual spectra resulting in 
a total of 3 INT spectra. 

Twelve additional blue spectra were obtained on the 
nights of 7-8 and 10-11 December 2008 with the Interme- 
diate dispersion Spectrograph and Imaging System (ISIS) 
double-arm spectrograph on the William Herschel Telescope 
(WHT). Here we employed the 600B and 1200B gratings on 
the blue arm resulting in wavelength coverages 3800-5400 
A and 3700-4500 A respectively. A 1" slit was selected yield- 
ing spectral resolutions in the range 65-115 km s~ . Mean- 
while, the 600R and 1200R gratings were mounted on the 
ISIS red arm, centered at different wavelengths, to yield si- 
multaneous spectra in the Ha and Call NIR triplet. Nightly 
averages were produced for the different spectral configura- 
tions. 

The High Efficiency and Resolution Mercator Echelle 
Spectrograph (HERMES) on the 1.2m MERCATOR (MT) 
telescope was also used to obtain three nightly 1800 sec- 
ond spectra on the nights of 30-31 October, 1-4 and 6-8 
November 2009. HERMES is a fiber-fed spectrograph and 
we employed the high resolution mode which yield a resolv- 
ing power R = 85000 accross the entire optical range be- 
tween 3770-9000 A. The signal-to-noise in the HERMES 
spectra is significantly lower than in the INT and WHT 
data so we decided to co-add the three individual spectra 
obtained every night resulting in a total of 9 Mercator spec- 
tra. 

Twenty two high-resolution echelle spectra were also ob- 
tained in the period Jan 20-Apr 8 2010 using the fiber-fed 
STELLA Echelle Spectrograph (SES) of the 1.2m robotic 
STELLA-I (ST) telescope at the Observatorio del Teide 
in Tenerife. The spectra cover the wavelength range 3870- 
8800 A with increasing inter-order gaps starting at 7200 
A. The spectrograph provides an effective resolving power 
R =55000. One spectrum was obtained per night and the 
exposure time was set to 1800 s. 

Finally, MWC 148 was also observed with FRODOspec 
at the LT between 4 Sep 2010 and 3 May 2011 using the 
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Table 1. Log of the MWC 656 observations. 



Date 






Telescope 


Spcct. 


Rangc 


Number of 


Exp. Time 


Dispersion 














spectra 


(seconds) 


(A pix- 1 ) 


23 Apr 


- 28 


Jul 2011 


LT 


3900- 


-5215 


32 


600 


0.35 


23 Apr 


- 28 


Jul 2011 


LT 


5900- 


-8000 


64 


290 


0.80 



Table 2. Log of the MWC 148 observations. 



Date 


Telescope 


Spcct. Range 


Number of 


Exp. Time 


Dispersion 






A 


spectra 


(seconds) 


(A pix" 1 ) 


20 Oct 2008 


INT 


3900-5500 


10 


90 


0.63 


21 Oct 2008 


INT 


3900-5500 


5 


90 


0.63 


22 Oct 2008 


INT 


3900-5500 


10,10 


120,90 


0.63 


07 Dec 2008 


WHT 


3800-5400 


3 


20 


0.88 


07 Dec 2008 


WHT 


8300-8960 


3 


20 


0.99 


08 Dec 2008 


WHT 


3800-5400 


3 


20 


0.88 


08 Dec 2008 


WHT 


5555-7070 


3 


20 


0.99 


10 Dec 2008 


WHT 


3700-4500 


3 


25 


0.45 


10 Dec 2008 


WHT 


8335-9075 


3 


15 


0.48 


11 Dec 2008 


WHT 


3700-4500 


3 


25 


0.45 


11 Dec 2008 


WHT 


8335-9075 


3 


15 


0.48 


30 Oct-8 Nov 2009 


MT 


3770-9000 


27 


1800 


0.02 


20 Jan-8 Apr 2010 


ST 


3870-8800 


22 


1800 


0.10 


4 Sep 2010-3 May 2011 


LT 


3900-5215 


40 


600 


0.35 


4 Sep 2010-20 Oct 2010 


LT 


5900-8000 


48 


193 


0.80 


24 Oct 2010-3 May 2011 


LT 


5900-8000 


47 


200 


0.80 



same configuration as for MWC 656 (see above). One 600s 
spectrum and two 290 s or three 193 s red spectra were 
obtained per night. A full log of the observations is presented 
in Table [2 

The INT, WHT and MT spectra were reduced us- 
ing standard techniques including debiasing and flatfield- 
ing. The spectra were subsequently extracted using opti- 
mal extraction techniques i n order to o ptimize the signal- 
to-noise ratio of the output IjHorne ' 1986). Frequent observa- 
tions of comparison arc lamp or hollow cathode lamp images 
were performed in the course of each run and the pixel-to- 
wavelength scale was derived through polynomial fits to a 
large number of identified reference lines. The final rms scat- 
ter of the fit was always <l/30 of the spectral dispersion. 
The automatic pipeline products were used for the ST and 
the LT data. 

The spectral type standard HR 2479 (BO III) was also 
observed with all different telescopes and instrument config- 
urations for the purpose of computing the rotational broad- 
ening of the companion star. It was selected bec ause of its 
low p rojected rotational velocity of 50 km s -1 l|Abt et al.l 



3 DESCRIPTION OF THE SPECTRA AND 
ROTATIONAL BROADENING 

All the spectra were rectified by dividing a low order spline 
fit to the continuum. For comparison, Figure [T] displays the 
averaged normalised LT spectra of MWC 148 and MWC 656. 
Both targets show similar spectra, with very broad and shal- 
low photospheric absorption lines. Strong Balmer emission 
is seen out to H/3 in MWC 656 and to H 7 in MWC 148. 
The H/3 line shows a double-peaked profile characteristic of 
circumstellar discs. Several Fe II disc emission lines are also 
detected. In summary, both spectra are typical of early-type 
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Figure 1. Average LT spectra of MWC 148 and MWC 656 in 
the blue (top panel) and red (bottom panel) spectral regions. The 
bottom spectra show the B0 III star HR 2479 for comparison, 
which has an intrinsic broadening of 50 km s — 1 . 
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Be stars with well develo ped discs fed by stable mass loss 
l|Porter fe Riviniusl 2003)) . Deta iled spectral analys i s per - 
formed by IWilliams et"aTl l|2010l ) and lAraeona et alJ <|2010h 
support a spectral classification B3 IVne+sh for MWC 656 
and BO Vpe for MWC 148. In the remaining of this paper 
we have adopted the spectral classification and stellar pa- 
rameters derived in these works. 

Projected rotational velocities of 430 km s -1 
(iGutierrez-Soto et alJl2007T ) and 500 km s" 1 (|Aragona et all 
l201Ch h ave been repo rted for MWC 148 whereas 262±26 
km s " 1 jYudin ll20oih and 300±50 km s _1 jWilliams et all 
bold ) for MWC 656. We decided to estimate the rotational 
broadening usini fr om our own spect ra following the 
technique outlined in iMarsh et al.l (1 19941 ). which basically 
subtracts broadened templates from the average spectrum 
of the target and searchs for the lowest residual. Therefore, 
we first rebinned our LT spectra of the two stars and the 
template HR 2479 into a uniform velocity scale of 23 km 
s _1 per pixel. The template was subsequently broadened 
from 100 to 500 km s -1 in steps on 5 km s" 1 using a 
Gray rotational profile (|Grav|[r992i ) with a limb darkening 
coefficient e = 0.33. The broadened versions of the template 
were multiplied by a scaling factor < 1 to account for 
dillution due to extra sources of continuum light (such as 
the equatorial Be disc) and subtracted from the averaged 
spectra of the two targets. The subtraction was restricted 
to the blue spectral range 4400-5200 A, where the photo- 
spheric He I, He n and metallic lines are most prominent. 
The IS, Balmer and Fe n emission lines were all masked in 
the process. A % 2 test on the residuals of the subtraction 
yields optimum rotational broadening of 370 km s _1 for 
MWC 148 and 342 km s _1 for MWC 656, with a formal 
uncertainty of ± 10 km s _1 . The rotational velocity of the 
template star needs to be added quadratically to these, 
resulting in ^sini = 373 km s" 1 for MWC 148 and 346 
km s _1 for MWC 656. The same analysis was repeated for 
MWC 148 using the other instrument configurations and 
we always obtain values in the range 355-381 km s _1 . 



4 ANALYSIS OF MWC 656 

Be discs in X-ray binaries typically display long-term super- 
orbital variability associated to activity episodes, but some- 
times aIso_varjabJIitvm£du with the binary orbital pe- 
riod l|Zamanov et al. 1 1 19991 ). The strong Ha emision is the 
best tracer of disc variabiliy and hence we decided to anal- 
yse several line parameters in MWC 656 and test whether 
they are modulated with the claimed 60.37 d orbital pe- 
riod ^Williams et al.ll2010l ). First, we measured the equiva- 
lent width (EW hereafter) by integrating the Ha emission 
profile in the 64 individual spectra. In addition, the FW HM 
and velocity shift of the line centroid were extracted through 
simple Gaussian fits. The left panel of Figure [5] presents 
the evolution of the Ha parameters with time. Because our 
base line extends over just 96 days we cannot probe for the 
presence of a 60.37 d modulation through a period search 
analysis. 

The EW is fairly stable during our observations with a 
mean of ~23 A and a smooth increase toward the end of our 
observing window. A peak is detected at HJD 2455725 but, 
unfortunately, we lack the necessary time coverage to check 




0.4 0.6 
Orbital phase 

Figure 3. Radial velocity curve of MWC 656 folded on the 60.37 
d period. The best-fitting solution, using an eccentric orbit with 
e = 0.4, is overplotted. The dotted point (at phase ~0.0) has been 
masked from the fit. Phase is set to HJD 2453243.3. 



for repeatability with a 60 d period. On the other hand, the 
FWHM and centroid velocity do display significant vari- 
ability, with two maxima elapsed by ~60 d, in good agree- 
ment with the photometric modulation of I Williams et al.l 
(2010). The behaviour of the Ha line parameters with the 
60.37 d period is presented in the right panel of Figure [5] 

Next, we attempted to measure radial velocities from 
the photospheric lines of the Be star. This is complicated by 
the fact that these lines are extremely broad and, in most 
cases, blended with disc emission lines e.g. the He I lines at 
4922 A and 5015 A partly overlap with Fe n emissions at 
4924 A and 5018 A (see Figure [TJ. After careful comparison 
of the average spectrum of MWC 656 with the broadened 
template HR 2479, we decided to restrict the radial velocity 
analysis to three spectral windows covering the He I lines 
4471 A, 4713 A and 5048 A. We measured radial veloci- 
ties by cross-correlating each spectrum of MWC 656 with 
a template formed from the average of the entire database. 
We find that this yields better results than using the broad- 
ened version of HR 2479 as template. The radial velocity of 
the template was determined by fitting a Gaussian to the 
core of the He I 4471 A line and it was added to the veloc- 
ities obtained from the cross-correlation. Figure [3] presents 
the final radial veloc ities folded on the 60.37 d period of 
IWilliams et~aH <|2010l ). The plot reveals a sine-like modula- 
tion which hints for a moderate eccentricity. 

The radial velocity curve was subsequently modelled 
with an eccentric orbital solut ion using t he Spectroscopic 
Binary Orbit Program (SBOP, lEtze1ll985D . The orbital pe- 
riod was fixed to 60.37 d. and phase arbitrarily set to 
HJD 2453243.3 which corresponds to the epo ch of maximum 
brigh tness in the photometric light curve ( Wi lliams et al.l 
2010) . Individual points were weighted proportionally to 
1/cr 2 , where a is the radial velocity uncertainty. In the solu- 
tions we adjusted the following orbital parameters: argument 
of the periastron (ui), systemic velocity (7), phase of the pe- 
riastron (</> pe ri) and velocity semiamplitude (K a pt)- The so- 
lution does not converge if the orbit eccentricity (e) is left 
free. Therefore, we attempted several fits fixing the eccen- 
tricity by hand and find a minimum rms for e = 0.4. Fitting 
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Figure 2. Left: time evolution of the Ha line parameters in MWC 656. Both the FWHM and centroid velocities show significant 
variability with maxima roughly separated by ~60 d. Right: Same but folded on the 60.37 d period of Williams et al. (2010). Phase 
has been set to the time of maximum optical brightness at HJD 2453243.3. The vertical dotted line denotes the phase of the periastron 
(see text). 



Table 3. Orbital solution for MWC 656 and MWC 148. 



Parameter 


MWC 656 


MWC 148 


P orb (days) 


60.37 (fixed) 


321 (fixed) 


T (HJD-2,450,000) 


3243.3 (fixed) 


4857.5 (fixed) 


e 


0.4 (fixed) 


0.83±0.08 


w (deg) 


71±23 


129±17 


7 (km s — - 1 ) 


-2.8±9.4 


48.3±8.9 


^pcri 


0.74±0.05 


0.967±0.008 


K opt (km s _1 ) 


41.7±6.8 


22.0±5.7 


ai sini (Rq) 


45.6±7.3 


77.6±25.9 


f(M) (M ) 


S5+ - 20 


06 +0 15 

u - uo -0.05 


cr (km s _1 ) 


20.3 


12.8 



tests using several orbital configurations indicate that an 
uncertainty of ±0.1 for e = 0.4 is appropriate. Note that we 
have masked the one discrepant point at phase from the 
fit. A close look to this spectrum shows a flatter He I 4471 
A profile than the rest, possibly due to disc emission. In any 
case, this point has a marginal impact in the final solution 
since it slightly reduces the eccentricity to 0.35 while the re- 
maining orbital parameters are virtually unaffected. Table[3] 
presents our final best-fitting parameters. 



5 ANALYSIS OF MWC 148 



In this section we repeat the previous analysis for the star 
MWC 148. The left panel in Figure g] presents the be- 
haviour of the main Ha parameters over our baseline of 
nearly three years. The LT data reveal a smooth sinusoidal 
modulation in all three parameters with a timescale which, 
at first glance, seem s consi stent with the 321 d X-ray period 
of iBongiorno et "all (|201ll ). This is better depicted in the 
right panel of Figure SI where the Ha parameters are folded 
on the 321 d period. The plot also shows that the EWs of 
the MT and, in particular, the ST spe ctra are systematically 
lower than the rest. For comparison, lAraeona et al.l (|2010l ) 
report £W=52.3 A on JD 2454757-92. The FWHM and 
centroid velocities are again consistently lower in the MT 
and ST spectra. A plot of the average Ha profiles demon- 
strates that the line becomes weaker and narrower in the 
MT and ST spectra because of the gradual appearence of 
a broad absorption throat with a full-width-zero-intensity 
FWZI ~2000 km s" 1 (see bottom panel in Figure HJ. This 
is characteristic of the faint states of Be stars which are 
likely related to episodes o f reduced circumstellar envelope 
(|Grundstrom fe Gies 1120061 ). Broad absorptions are also de- 
tected in H/3 and the myriad of Fe n lines which plague the 
blue spectral range such as Fe n 4924 A and 5018 A (see top 
panel in Figure [5]). These absorptions will likely contami- 
nate most of the photospheric lines. Therefore, we decided 
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Figure 4. Left: time evolution of the Ha line parameters in MWC 148. Different symbols indicate different datasets as follows: asterisk 
(WHT), open triangles (MT) open circles (ST) and filled circles (LT). Note the sinusoidal modulation in the LT data with a timescale 
of «300 d. Errorbars are not plotted because they a re always smaller than the symbol size. Right: same as in the left panel but folded 
on the 321 d X-ray period of lBongiorno et al.l j201ll ). Phase has been set to HJD 2454857.5. Periastron is marked by a vertical dotted 
line (see text). 



to exclude these spectra from the study of the radial velocity 
curve of the Be star. 

As in MWC 656, radial velocities were obtained by 
cross-correlating every blue INT/WHT/LT spectrum with a 
template produced from the average of all. Prior to this, each 
spectrum was rebinned into a uniform velocity scale of 23 km 
s _1 per pixel. Cross-correlation was performed over spectral 
windows covering the He I lines 4471 A, 4713 A and 5048 
A. Additional regions containing lines of Si ill (4568 & 4575 
A), C m/0 II (4639-4650 A) and He II 4686 A were also in- 
cluded. These absorption features are stronger in MWC 148 
than in MWC 656 because of its earlier spectral type. The 
rest velocity of the template was again determined through 
a Gaussian fit to the core of the He I 4471 A line and it was 
subsequently added to the cross-correlation velocities. The 
radial velocity curve, folded on the 321 d period is displayed 
in Figure [6] In spite of the limited phase coverage and ve- 
locity scatter, a sharp velocity minimum is detected. The 
narrowness of this feature prompts for a large eccentricity. 
We have attempted to model this radial velocity curve with 
the SBO P fitting code, fix i ng th e orbital period to 321 d. 
Following iBongiorno et aD ()201ll ). we have defined phase 
as HJD 2454857.5 which was arbitrarily set to the date of 
the first Swift/XRT observation in their paper. The best- 
fitting orbital elements are presented in the second column 
of Table O 



6 DISCUSSION 



In this paper we have presented an extensive spectroscopic 
database of MWC 656 and MWC 148, the optical counter- 
parts of the candidate 7-ray binaries ACL J2241+4454 and 
HESS J0632+057 respectively. Our data show evidence for 
long-term modulation in the main Ha parameters, consis- 
tent with the proposed orbital periods of 60.37 d and 321 d 
respectively. These modulations are likely produced by tidal 
instabilities in the circumstellar envelopes triggered by the 
motion of the compact star in an eccentric orbit. Circum- 
stantial support for the orbital nature of the Ha variabil- 
ity is given by the EW values. The EW of the Ha line in 
Be/X-ray binaries provides a good estimate of the size of 
the circumstellar disc and, because this is truncated by the 
tidal torques of the co mpact star, a simple relation iJWXHa) 
oc P 4/3 is expected (jReig I I201J ). It should be noted that 
this correlation holds for the maximum EW observed over 
a long length of time. We measure EW = 25 A for MWC 656 
and 56 A for MWC 148 which, accor ding to the empirical 
P orb - EW{Ra) diagram in Fig. 15 of iReig I (|201lf ). suggest 
Porb ~ 90 and 250 d respectively. Given the scatter of the 
figure, these are fully consistent with the orbital periods pro- 
posed in the literature. 
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Figure 5. Average spectra of MWC 148 showing the evolution 
of several circumstellar disc profiles. The emission lines weaken 
because of the appearence of broad absorption throats, especially 
clear in the ST spectra. 



6.1 Masses of the compact stars 

Further, we have detected radial velocity variations in the 
photospheric lines of both MWC 656 and MWC 148. As- 
suming that their orbital periods are 60.37 d and 321 d, we 
produce phase folded radial velocity curves which allows us 
to constrain their orbital elements for the first time. In par- 
ticular, we determine the mass function of the compact ob- 
ject to be f(M) = M c 3 sin 3 ij (M c + M opt ) 2 = 0.35±g"?§ M 
for MWC 656 and 0.06± ; | M Q for MWC 14. In this equa- 
tion, M c and M op t stand for the masses of the compact and 
the optical star respectively. Unfortunately, no pulsations 
have been detected from the compact stars yet and hence we 
cannot provide a full solution to the stellar masses and the 
nature of the unvisible companions at this point. However, 
some speculations can be made based on our determination 
of the rotational broadenings usini. 

We measure v sin i ~ 346 km s- 1 for MWC 656 and 
~ 373 km s _1 for MWC 148. A lower limit to the incli- 
nation can be derived from the condition that the optical 
companion should not exceed 0.9 times th e critical rota - 
tional velocity w C rit. According to Table 2 in lYudin I l|200lf ). 
«crit ~ 565 km s _1 for a B0 V which implies i>47° for 
MWC 148. On the other hand, v clit ~ 420 km s _1 for a 
B3 IV and hence i>66° for MWC 656. These estimates 
obviously assume that the Be star's spin axis and the or- 
bital axis are aligned, which may not be true considering a 
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Figure 6. Radial velocity curve of MWC 148 folded on the 321 d 
period. The best eccentric orbital solution is overplotted. Phase 
is set to HJD 2454857.5. 



possible kick during the supernova explosion that forms the 
compact object. 

In addition, the presence of double peaked emission pro- 
files in our spectra also hints for a moderately high incli- 
nation. In particular, the Fe n lines are optically thin and 
their profiles reflect t he Keplerian rota tion in the emitting 
part of the Be disc <|Hanuschiklll996h . Therefore, we can 
use their FWZI to estimate the binary inclination through 
FWZI/2 x sini = (GM opt /R opt ) 1/2 . We measure FWZI 
(A5018) ~1300 km s" 1 in MWC 148 and, using M opt ^13.2 - 
19.0 M and R opt ~6. 0-9.6 R Q from lAragona et all i|2010h . 
we find i ~71-90° . Regarding MWC 656 we measure FWZI 
(A5018) ~1000 km s" 1 which, com bined with M op t -5.8 - 
9.8 M Q and R opt ~4.7-8.5 R from I Williams et all < |2010h . 
yield i ^70-76°. These estimates, however, should be re- 
garded as mere upper limits because the region where the 
Fe II lines are produced is unlikely to extend down to the 
surface of the Be star. 

We also note that the Fe II profiles do not exhibit 
deep central absorptions below the continuum, character- 
istic of shell stars. Shell lines are thought to be produced 
by obscuration of the star b y the circumstel lar disc and 
are only observed at i > 80° (Hanuschik 1996), which sets 
a rough upper limit to the inclination in both MWC 148 
and MWC 656. In summary, crude values of i ~47-80° for 
MWC 148 and i ~67-80° for MWC 656 seem the most plau- 
sible given the available constraints at this point. Bringing 
the c onstraints to th e inclin atio n and the Be st e llar m asses 
from IWilliams et al.l i|2010l ) and lAragona etafl (|2010T ) into 
the mass function equation results in M c ~2.7-5.5 M for 
MWC 656 and 1.3-7.1 M for MWC 148, although these 
numbers should be treated with caution because of the large 
uncertainties involved in this calculation. First, the binary 
inclinations are loosly constrained and need to be refined. 
Further, our orbital solution is based on a limited data set 
which requires confirmation through a more extended base- 
line. In particular, observations of MWC 148 around the pe- 
riastron at phase are strongly encouraged. Unfortunately, 
the long orbital period coupled to the annual solar cycle 
means that the first opportunity for ground observations 
will not happen until March 2015. 

The wide range of masses allowed by our orbital solu- 
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tions prevent us from deciding whether the compact objects 
in MWC 148 and MWC 656 are neutron stars or black holes. 
However, their long orbital periods place them among the 
widest 7-ray binaries, only after PSR B1259-63, and this 
offers an opportunity to probe the nature of the compact 
star through the detection of radio pulsations. Previous at- 
tempts in LS 5039 and LS I +61 303 have failed presum- 
ably due to strong free-free absorption by the dense stellar 
winds (see McSwain et al. 2011 and references therein). In 
principle, the detection of radio pulses might be possible 
in both MWC 148 and MWC 656 around apastron phases 
due to their wider orbits. X-ray pulses are also expected 
to arise through synchrotron emission in the magnetosphere 
of a rotationally-powered pulsar or from an accreting neu- 
tron star. Deep searches for X-ray pulsations have been per- 
formed in both LS 5039 and LS I +61 303 with null results 
(See Rea et al 2010, 2011 and included references). This 
implies than the putative pulsars are spinning faster than 
~5.6ms, the pulsar beam is pointing away from our line of 
sight or X-ray pulsations are r estricted to a lim ited range 
of orbital phases. We note that lRea et al. 1 12010) could not 
either find pulsations in MWC 148 using Chandra data dur- 
ing one outburst nor in XMM — Newton data during qui- 
escence. These observations had a time resolution of only 
199.2 ms and, according to our ephemeris, were performed 
at orbital phases 0.33 and 0.46, respectively. Therefore, the 
XMM ~ Newton data coincide with the apastron passage, 
where Compton scattering is lower and hence more favor- 
able for the detection of X-ray pulses. More higher time 
resolution observations of both MWC 656 and MWC 148 
at different orbital phases are urgently required to test the 
nature of their compact stars. 

6.2 Orbital variability 

Aside from the stellar masses, our orbital solutions do pro- 
vide interesting information on the geometry of the 7-ray 
binaries and properties of the observed X-ray /VHE radia- 
tion. Figure [7] presents the relative motion of the compact 
star around the optical companion for both MWC 656 and 
MWC 148 as seen from above i.e. i = 0°. The figure was pro- 
duced using the stellar mas s es and radii from I Williams et al.l 
(|2010h and lAragona et ail (|2010h and assuming a 1.4 M 
compact star. Our ephemeris imply that the X-ray outbursts 
in MWC 148 are delayed by ~100 d or 0.3 orbital phases 
with respect to the periastron passage. This is remarkably 
similar to LS I +61 303 where X-ray outbursts are observed 
~0.2 orbital phases passed the clos est approach of the stellar 
components (jCasares et al.ll2005ah . 

Regarding MWC 656, the maximum of the optical mod- 
ulation takes place ~0.25 orbits after periastron. A sim- 
ilar phase shift of the photometric maximum is seen in 
LS I +61 303 (|Mendelson fc Mazehl 1 19891 ; iParedes et ail 
Il994h . In addition, the right frame in Figure [5] shows that 
the EW of the Ha line peaks ~0.3 orbital phases passed 
periastron. This behaviour seems consistent with an sce- 
nario where the modulation in both the continuum flux 
and the (opticaly thick) Ha flux is caused by the visibil- 
ity of a one-arm tidal wave in the circumstellar disc, trig- 
gered by the c l ose pa ssage of the neutron star. For instance, 
Okaz aki et al.l l|2002h present SPH simulations of a Be/X- 
ray binary with P or b=24.3 d and e = 0.34 which is not 



too far from the parameters that we derive for MWC 656. 
Their simulations (see Figures 10 and 11 in the paper) show 
that a spiral density wave, excited by the periastron pas- 
sage, is fully developed aroun d phas e ~0.25 past periastron. 
Furthermore, IWilliams et alj (|2010h observed a rapid vari- 
ability of the V/R ratio of the H7 emission and suggested 
it could be associated with the changing tidal effects of the 
companion on the disk, which is specially strong near perias- 
tron. Indeed, these observations correspond to orbital phases 
0.82 + 0.02 and 0.84 + 0.02 which according to our ephemeris 
are just ~0.08 and ~0.10 in phase after periastron. At phase 
~0.8 we also observe a rapid transition of the Ha centroid 
from negative to positive velocities, coincident with a peak 
in the FWHM (see right frame in Figure [2} ■ 

In view of the new orbital ephemeris of MWC 656 it 
is also interesting to put in context the multiwavelength 
observations of the source. First of all, the AGILE detec- 
tion between 25 July 2010 at 01:00 UT and 26 July 2010 
at 23:30 UT l|Lucarelli et al. I I2010T ). corresponding to an 
orbital phase range of 0.77-0.80, took place just after (or 
even during) periastron. A maximum of the GeV emis- 



sion takes place during periastron in LS 5039 (Abdo et al 



2009a) and soon after it in LS I +61 303 (|Abdo et al 



2009b). The fact that the AGILE Rare in MWC 656 happens 



very close to the inferior conjunction of the compact object 
(w = 71 ± 23 deg), coupled with its high orbital inclina- 
tion, may explain the weak GeV emission through a reduc- 
tion in the Compton scat tering cross-section for large angles 
l|Khangulvan et al. 1 120081 ). All these facts reinforce the asso- 
ciation between the Be star MWC 656 and the gamma-ray 
source AGL J2241+4454 and, even with the lack of 7-ray 
detections aside from the flare seen by A GILE during a sin- 
gle periastron passage, give strong support to the idea that 
MWC 656 is a new 7-ray binary. 

On the other hand, MWC 148 has only shown signifi- 
cant VHE radiation simultaneous to the X-ray flare episodes 
which occur ~0.3 phases past periastron. Ind eed, according 
to ou r ephemeris, the VERITAS upper limits (|Acciari et al.l 
2009) were obtained during the phase intervals 0.6-0.72, 
0.86-0.93 and 0-0.0 1, corresponding to low levels of quies- 
cent X-ray emission (|Bongiorno et al]|201ll ). Conversely, the 
H.E.S.S. discovery observations occured at phases 0.32 and 
0.36-0.46, coincident with periods of X-ray flares. Further 
detections of MWC 14 8 at VHE during X-ra y active phases 
have b een reported bv lOng et al.1 (|2011 ) and lMariotti et all 
|201ll ). 



6.3 The class of 7-ray binaries 

The family of 7-ray binaries is growing fast with the re- 
cent discovery of AGL J2241+4454, HESS J0632+057.1 and 
1FGL J1018. 6-5856. These add up to the group of "classic" 
7-ray binaries PSR B1259-63, LS 5039 and LS I +61 303. 
Table 0] summarizes some of their main physical parameters. 
They all contain late O-Be stars with strong winds and stel- 
lar radii in the range ~7-10 Rq. The table shows a possible 
correlation between the eccentricity and the orbital period 
and this is illustrated in figure [8] A least-squared linear fit 
yields: 



e = 0.206(19) x logP orb + 0.233(60) 



(1) 
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Figure 7. Left:relative orbit of the compact object around the optical star in MWC 656 as seen from above the orbital plane. Phases 
of periastron and apastron are indicated and joined by a dashed line which depicts the major axis of the orbit. Stars mark 0.1-phase 
intervals from periastron whereas open triangles indicate the phases of the inferior conjunction and superior conjunction of the compact 
star. The coordinates are in units of the orbital semimajor axis. Right: same for MWC 148. Note that the semimajor axis is 0.64 AU in 
MWC 656 and 2.38 AU in MWC 148 



For comparison, Figure [8] also presents 35 high-mass X-ray 
bin aries (HMXBs ) with X-ray pulsars from lLiu et al. I (|2006T ) 
and lMartin et all (|2009l ). HMXBs tend to lie below the linear 
regresion of 7-ray binaries, i.e. for a given orbital period 
they have lower eccentricities. The only exceptions are 4U 
1850-03, 0535-668 and the SMC binary J0045-7319, all 
with e > 0.8. In particular, the latter two have P sp m < 
1 s and hence a high spin-down power which makes them 
potential candidates for 7-ray emission. 

The reason behind the P or b — e correlation in 7-ray bi- 
naries stems from the fact that small stellar separations are 
required to trigger VHE radiation and this can only hap- 
pen for long orbital periods if the eccentricity is large. Con- 
versely, for short orbital periods high eccentricities would 
imply that the compact object passes through the stellar 
atmosphere of the companion, clearly an unstable situation. 
The separation of the two stars at periastron a(l — e) and 
apastron a(l + e) is also listed in Tabled! assuming that 
the compact star is a 1.4 Mq neutron star. As expected, the 
most compact binaries LS 5039, 1FGL J1018.6— 5856 and 
LS I +61 303 show persistent VHE and/or HE radiation 
modulated with their orbital periods whereas only transient 
radiation is detected in PSR B1259— 63 around periastron. 
This seems to imply that stellar separations < 1 a.u. or < 25 
i? opt are required for the efficient production and modula- 
tion of the HE/ VHE radiation but the picture is certaintly 
more complex. For instance, MWC 656 has been detected 
up to GeV energies only once during a flare episode at peri- 
astron whereas MWC 148 is regularly detected at TeV ~0.3 
orbital phases past periastron. Aside from sensitivity limita- 
tions of different satellite/telescopes, other processes such as 
photon-photon annihilation and the subsequent production 
of electron-positron pairs are likely responsible for the wide 
phenomenology observed. In particular, the binary inclina- 




Figure 8. Eccentricity versus orbital period for 7-ray binaries 
(filled circles) and high mass X-ray binaries from the catalogue of 
Liu et al. 2007 (open circles). A linear fit to the 7-ray binaries is 
overplotted. 



tion and argument of periastron are important parameters 
for the attenuation of the VHE rad iation and modelling the 
observed amplitudes (Dubus 2006). 

It is well established that decretion discs in Be/X-ray 
binar ies are truncated by the tida l torque of the neutron 
star (jNegueruela fc OkazakT1l200ll ). Since the EW of the 
Ha emission is a good proxy of the disc size this results in 
the well known corr elation between P or b and the maximum 
observed EW(Ha) (|Reig et al. Ill997l ; iReig Il201ll ). Most 7- 
ray binaries contain Be primaries and it is interesting to 
test whether a similar correlation holds. Table [4] lists the 
range of EW(Ha) values reported in literature for the 7- 
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Figure 9. Maximum EW of the Ha emission line versus orbital 
period for 7-ray binaries (filled circles) and Be/X-ray binaries 
(open circles). The latter include Milky Way and SMC binaries 
from Reig (20f 1). A linear fit to the 7-ray binaries is overplotted. 
Note that LS 5039 is not included because the primary is not a 
Be star. 

ray binaries. Significant long-term and orbital variability is 
typically seen in Be/X-ray binaries and the same is observed 
in the 7-ray binaries. In this paper we detect a maximum 
EW for MWC 656 at ~ 0.3 phases past peri astron, in agree- 
ment with what is observed in LS I +61 303 (|McSwain et al.l 
l2010h . On the other hand, the maximum EW for MWC 148 
occurs at periastron. In any case, only the maximum EW 
values are driven by the disc truncation radius. Figure [9] 
plots the maximum EW(Ha) measured in 7-ray binaries 
with Be primaries versus P or b and we note that they follow 
the linear regresion: 

EW = 37.8(5.0) x log P orb - 38.3(11.4) (2) 

For compar ison, Figure EJ] also shows EWs of Be/X-ray bi- 
naries from iReigl (|201 lh . As in the case of the eccentricity, 
we observe that 7-ray binaries define the higher envelope of 
the distribution i.e. for a given orbital period Be/X-ray bina- 
ries tend to have smaller decretion discs than 7-ray binaries. 
This is somehow surprising in a millisecond pulsar scenario 
because it implies that the relativistic pulsar wind does not 
"erode" significantly the circumstellar disc even during peri- 
astron passage. However, it is also true that, given the large 
eccentricities of 7-ray binaries, they only spend a short frac- 
tion of time close to the Be star, perhaps preventing efficient 
truncation. Interestingly, equations 1 and 2 provide a new 
tool for a rough estimate of two fundamental parameters 
(such as P or b and e) in Be/7-ray binaries simply from the 
EW of the Ha emission. This can be tested with future 
observations of newly discovered 7-ray binaries. 



7 SUMMARY 

We have reported optical spectroscopy of MWC 656 and 
MWC 148, the optical counterparts of the 7-ray sources 
ACL J2241+4454 and HESS J0632+057 respectively. Our 
data show that the main Ha parameters are modulated with 
the 60.37 d optical photometric period in MWC 656 and 



with the 321 d X-ray period in MWC 148. This is likely 
produced by the visibility of a tidal wave in the Be disc, 
triggered by the close passage of a compact star. In addi- 
tion, we present radial velocity curves for the two stars and 
the first constraints to the binary parameters. Both stars 
display similarities with LS I +61 303, where photometric 
maxima and X-ray outburst are delayed by ~0.3 phases past 
periastron. When compared to other 7-ray binaries, we find 
that the eccentricity is correlated with the orbital period. 
This is explained by the small stellar separations needed for 
triggering the 7-ray activity. The maximum EW of the Ha 
emission is also correlated with orbital period. This is rem- 
iniscent of Be/X-ray binaries and suggest that the size of 
circumstellar discs in 7-ray binaries is tidally truncated by 
the compact companion. 
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